% 创建时间：2025-06-22 00:13:59
% 作者：兪馳
clc;
clear;

H=diag([1,1,1],1),
for  i=2:4,
    H^i,
end